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SUMMARY 


Measurements of streamwise velocity fluctuation and surface pressure 
fluctuation spectra and wavespeeds are reported for a well -documented 
separating turbulent boundary layer. A two-parallel -wire hot-wire anemom- 
eter probe was used to measure the velocity fluctuations while two sensi- 
tive instrumentation microphones were used to measure pressure fluctuations 
through pinhole apertures in the flow surface. Because a portion of the 
acoustic pressure fluctuations are the same across the nominally two- 
dimensional turbulent flow, it is possible to decompose the two microphone 
signals and obtain the turbulent flow contributions to the surface pressure 
spectra. ! i : 

The velocity spectra and wavespeeds behave as expected with the spectrum 
function F(n)~ n" 1 near the wall upstream of detachment and F(n) - n" 5 ^ 3 
in the higher frequency inertial subrange. The wavespeed or celerity 
increases at low frequencies to the local mean velocity at high frequencies, 
as observed by Strickland and Simpson (1973) for a separating turbulent 
boundary layer. 

The rms surface pressure fluctuation p' increases monotonically 

through the adverse pressure gradient attached flow region and the detached 

~2 

flow zone. Apparently p is proportional to the ratio of streamwise length 
scale to length scales In other directions a. For non-equilibrium 
separating turbulent boundary layers, a is as much as 2.5, causing p 2 to 
be higher than equilibrium layers with lower values of a. 

The maximum turbulent shearing stress t m appears to be the proper 
stress on which to scale p 1 ; p'/t M from available data shows much less 


i 


I 

! 


variation than when p' is scaled on the wall shear stress. In the present 
measurements P'/t m increases to the detachment location and decreases 
downstream. This decrease is apparently due to the rapid movement of the 
pressure-fluctuation-producing motions away from the wall after the begin- 
ning of intermittent backflow. A correlation of the detached flow data is 
given. 

The spectra §(w) correlate well when normalized on the maximum 

shearing stress t^. At lower frequencies for the attached flow §(u>) ~ 

-3 

while |(w) - u at higher frequencies in the strong adverse pressure 
gradient region. After the beginning of intermittent backflow, §(<*) varies 

_3 

with w at low frequencies and w at high frequencies; farther downstream 

2 4 

the lower frequency range varies with u ' . 

The celerity of the surface pressure fluctuations for the attached 
flow increases with frequency as suggested by other measurements and 
agrees with the semi -logarithmic overlap equation of Panton and Linebarger. 
After the beginning of the separation process, the wavespeed decreases 
because of the oscillation of the instantaneous wavespeed direction. The 
streamwise coherence decreases drastically after the beginning of flow 
reversal . 
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NOMENCLATURE 


c 

d 

F(k),F(n) 

H 

K 

k 

L 

N 

n 

P.P 

P' 

q 

t 

R 

n 

S 

U,V,W 


u,v,w 

T ~2 T 
u ,v , w 


-uv 

U M ,U 

M T 


U ,U 
crt cn 


~s 

X,Y, Z 

AX , AY , A z 


speed of sound, mps 
diameter of microphone, m 

spectrum function of k or n; j£°°F(n)dn = f Q F(k)dk = 1 
6*/e, velocity profile shape parameter 
damping constant in equation (28) 

= 2 it/ X wavenumber, m”^ 

length scale in equation (25), m 

distance from wall to maximum shear stress location, m 
frequency, Hz 

mean and fluctuation wall pressures. Pa 
rms wall pressure fluctuation. Pa 
dynamic pressure, Pa 


position vector, m 

correlation coefficient for frequency n 
length of detached flow zone, m 
mean velocities in X,Y, Z directions, mps 
fluctuation velocities in X,Y, Z directions, mps 


p 

mean square velocity fluctuations, (mps) 

2 

kinematic Reynolds shear stress, (mps) 




, shear velocities based on maximum and wall 

P P 

shear stresses s mps 

to instantaneous and mean celerities or wavespeeds at 
k ’ k 


frequency n, mps 

Perry and Schofield velocity profile scale, mps 
streamwise, normal to wall, spanwise positions, m 
spacing of sensors in streamwise, normal to wall, spanwise 
directions, m 


iv 



Greek Symbols 


a ratio of streamwise length scale to length scales in other 

directions 

2 2 
y coherence function, |G^y| /G^Gyy 

Y square root of the coherence function 

A Perry and Schofield length scale, m 

6,6*,e boundary layer, displacement and momentum thicknesses, m 
k von Kdrmdn constant 

2 

v kinematic viscosity, m /S 

n Coles wake function parameter 

3 

p density, kg/m 

x shearing stress. Pa 

2 ^ 

$(u>) spectrum function, p' = f ^(aj)du) 

-CO 

phase 

a) 2irn, angular frequency, rad/S 


Subscripts 

i , j tensor directions 

n value associated with frequency n 

ref reference inlet value 

00 local freestream condition at outer edge of boundary layer 


v 



I. INTRODUCTION 


Noise generated by helicopter and turbomachine rotors is a nuisance 
that designers would like to predict and to minimize within other design 
constraints. Brooks and Schl inker (1983) reviewed some recent research 
on helicopter rotor noise and discussed the categories of noise sources, 
which include blade self-noise generation by strong adverse-pressure- 
gradient attached turbulent boundary layers and by separated turbulent 
boundary layers that accompany stall. 

Brooks and Hodgson (1981) showed that starting with given surface 

pressure fluctuation spectra and convective speeds, radiated noise due to 
♦ 

the turbulent boundary layer can! be predicted. Furthermore, if the 
surface pressure fluctuation spectra and convective speeds can be related 
to the turbulent flow structure, then turbulent boundary layer flowfield 
calculation methods can be used when designing rotors to estimate the 
needed surface pressure fluctuation information. 

Thus, a key requirement for this noise calculation procedure is 
knowledge relating the flowfield structure to the surface pressure 
fluctuation structure. Unfortunately there are few measurements of both 
flowfield structure and surface pressure fluctuation structure for a 
given flow, especially in the presence of adverse pressure gradients. 

Only recently have detailed flowfield measurements been made of a nomi- 
nally two-dimensional, steady freestream, separated turbulent boundary 
layer (Simpson et al . , 1981, a, b, c). This report presents information 
on the surface pressure fluctuations for this wel 1 -documented velocity 
flowfield. 
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After discussing some previous work on surface pressure fluctu- 
ations for zero and adverse-pressure-gradient and separating turbulent 
boundary layer flows in the next section, a review of the test wind 
tunnel and a summary of the nature of this separating turbulent boundary 
layer flowfield are given. In following sections the experimental 
instrumentation, techniques and results are discussed for the streamwise 
velocity spectra and celerities and the surface pressure fluctuation 
spectra and celerities. 

II. SOME PREVIOUS WORK ON SURFACE PRESSURE FLUCTUATIONS 

First, it should be stated at the beginning that static pressure 
fluctuations can only be measured non-intrusively at the wall in a 
turbulent boundary layer. Small pinhole openings in the surface have 
been used to measure fluctuations with microphones. Small piezoelectric 
transducers flush with the surface have also been used to avoid any 
influence of a pinhole on the flow. Sufficiently small sensing surfaces 
must be used to pick up small scale fluctuations. Flow disturbances are 
present in most wind tunnels, so measurements of the low frequency 
components of the wall pressure fluctuations are impossible to make with 
a single sensor. Willmarth (1975) reviewed many measurements, pointing 
out the shortcomings of each set of available data. Needless to say, 
there are significant differences between results obtained for approxi- i 
mately the same flow conditions by different investigators. 

As in all other practical turbulent flow problems, purely theoret- 
ical calculations of surface pressure fluctuations have not yet been 
done successfully. In incompressible turbulent flow the fluctuating 
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pressure p is related to the velocity fluctuations by the Poisson 
equation 


i!e = 

>4 


pq 


(D 


where the source term q is given in tensor notation by 


3U i 9u i a 2 

2 — - — + — 

3x • ax. 3x . 3x . 

J " \J 


(u .u . - U -U . ) 
1 J 1 J 


( 2 ) 


and and u^ are the mean and fluctuation velocity components in the x^ 
direction. The first term on the right side of this equation represents 
the turbulence-mean shear interaction while the second term represents 
the turbulence- turbulence interaction. For a wall-bounded flow, if 
contributions from surface integrals are neglected, then the fluctuating 
pressure at a point X on the wall is given by 


P ( X ) 


p 

2tt 


/ 

Y >0 


q(R s ) dV (R $ ) 
~ -- 

(X - R) 


(3) 


where the volume integration is at position over the entire half- 
space containing the flow. This equation indicates that while surface 
pressure fluctuations are produced from sources in a large region of the 
flow, contributions from various sources drop off rapidly with their 
distance from the point under consideration. Although several attempts 
have been made to calculate p from equation (3), Willmarth (1975) pointed 
out that it appears that such efforts suffer from the lack of accurate 
information about the fluctuating velocity field within the boundary 
layer. Consequently, the effects of adverse pressure gradients and 
separation on surface pressure fluctuations as suggested by calculations 
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are uncertain and need confirmation by experimental data. 

Panton and Linebarger (1974) calculated wall pressure spectra for 
zero pressure gradient and adverse pressure gradient equilibrium boundary 
layers that seem to describe the essential features observed from experi- 
ments. They used Coles' laws of the wall and wake for the mean velocity 
profiles. A scale-anisotropic model of the spatial correlation of v was 

I /O 

used together with the assumption that v was proportional to (-uv) ' . 

As in some earlier calculations, only the turbulence-mean shear interaction 
was modeled since the turbulence-turbulence interaction contributes only 
a very small portion to the mean-square value. 

Their spectral results showed larger contributions at higher Reynolds 
numbers for ks < 20. The contributions to the spectrum at these low 
frequencies are due to the outer region velocity and turbulence structure 
and depend upon the pressure gradient. An overlap region between the low 
frequency outer region contributions and the high-frequency near-wall 
viscous-dominated part of the spectrum varies with k~\ as shown by 
Bradshaw (1967) using dimensional analysis. For this overlap region, 
their calculation results can be approximated by 

^j^-xl.73 (a ) 0 ' 9 for jp < 0.06 (4) 

T w T 

Here a is the ratio of the streamwise length scale to length scales in 
other directions, which strongly influences the spatial correlation of v. 

At higher frequencies 
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^^-=0.0173 (jf-\ for {p > 0.1 (5) 

T w T T 

Note that these two equations are independent of Reynolds number and are 
scaled on the wall shearing stress. 

Because the low frequency part of the spectrum is Reynolds number 
dependent, the mean square pressure fluctuation increases with Reynolds 
number. The equation 


= 0.52 
T w 


0 9 V 

x Un / - 7 - / + 9.24) 


( 6 ) 


fits Panton and Linebarger's calculations for a zero pressure gradient 
with a = 1, 2 and 3, and a iColes' wake function parameter n = 0.6. In 
terms of the displacement thickness Reynolds number 


xi. U 


U S 

T 

V 


U 6 


- 65 ^ 


(7) 


1+n \ v 

(Coles and Hirst, 1969). 

The spectral measurements of Bull and Thomas (1976) showed that zero 
pressure gradient flow measurements with a pinhole microphone, such as 


Blake's (1970), are too high for frequencies higher than the high 

2 

frequency end of the overlap region, wv/U^ > 0.1. This effect extends 

9 

to wv/U > 2. The pinhole microphone interacts with the viscous flow 

T 

to amplify the frequencies in this range. The spurious amplification 

2 

appears to be a function of tov/U^ for the narrow range of test parameters 
that were used. Bull (1967), Lim (1971) and Schloemer (1967) used flush- 
surface piezoelectric pressure transducers that were relatively large 
so high frequency contributions to the mean square pressure fluctuation 
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were attenuated (Willmarth, 1975). Consequently, the rms pressure 
fluctuations were too low for these latter measurements, as shown in 
Figure 1. 

The spectral measurements of Bull and Thomas (1976) show a short 

12 o 

$(u>) ~ of region (0.05 < wv/U t < 0.2) and a higher frequency 4>(u>) - u'- 3 

2 

region (u>v/U x > 0.5). Their rms pressure fluctuation measurements 
with a piezoelectric transducer are still slightly too low because the 
sensing diameter was too large (44 < dll^/v < 65) to pick up the smallest 
scale turbulent motions (a U /v ~ 12). Thus the p */t values calculated 

Z T W 

by Panton and Linebarger for the zero-pressure-gradient case appear to 
be approximately correct. 

At the lower frequencies of the overlap region, the pressure fluctuation 
spectrum for the zero-pressure-gradient case scales on the outer variables 
and 6 *. The spectral level seems to be almost the same in various 
investigations with 

= -51 + 2 at w6*/U = 1 (8) 

— oo 

(Bull and Thomas, 1976; Burton, 1973; Willmarth, 1975; Schloemer, 1966; 

Panton and Linebarger, 1974). For 0.1 < ooSVU^ < 1, $(uj) levels off and 
is about constant in most investigations. For the lowest measureable fre- 
quencies, $(u>) generally increases slightly with w, although there is 
some scatter in data due to different low frequency disturbances in different 
test wind tunnels. 

For adverse-pressure-gradient equilibrium turbulent boundary layers, 
the low frequency part of the spectrum (ks < 20) scales approximately on 


10 log 1Q 


$ (u))U a 
q 2 6 * 


7 


El 

T wall 



1000 10000 100000 


Re 


6 * 


Figure 1: p'/fwall for some zero pressure gradient experiments. 

Solid lines: Panton and Linebarger, a=1.0 and tt=0.6 (1974); 
□: Bull and Thomas (1976); a: Lim (1971); o: Schloemer; 
Shaded region: Bull (1967) and Willmarth (1958). 
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the maximum shearing stress in the boundary layer (Mabey, 1982). This 
can be seen in Figure 10 of Panton and Linebarger (1974) for tt= 1.5, 3, 
and 6, together with the maximum shearing stress estimated from 

' fejf (1 - (9) 

2 2 

(East and Sawyer, 1979). Table 1 also shows that p /x^ from Panton and 

~2 ? 

Linebarger is nearly constant while p /x varies over an order of 

w 

magnitude. This scaling on x M is supported by equation (3) and the struc- 
ture of strong adverse pressure gradient boundary layers. 

Perry and Schofield (1973) proposed universal empirical correlations 
for the inner and outer mean velocity profile regions of strong adverse- 
pressure-gradient boundary layers near separation. Their work was based 
upon 145 mean velocity profiles taken from Coles and Hirst (1969), including 
both equilibrium and non-equilibrium profiles. The data of Simpson et al . 
(1977; 1981 a, b, c) upstream of any flow reversal agree with these cor- 
relations. The Reynolds shearing stress profile and the maximum shearing 
stress X|^ (which occurs in the middle of the boundary layer) play important 
roles in these correlations when x^ > 1.5 t . In other words, the same 
large-scale turbulent structures that produce the maximum shearing stress 
away from the wall are the structures that influence the mean velocity 
profile from near the wall to the outer edge. 

Perry and Schofield proposed that the outer region flow be described 
by 


U -U 

— — = f (%-) 

U s t 2V ’ 


A = 2.86 


with 



Table 1. Some pressure fluctuation results from adverse-pressure 
gradient turbulent boundary layers 
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"1 = u)5 */Uo 0 ; * denotes values estimated from equilibrium boundary layer equations for '-he laws of wall 
wake or equation (9); NA denotes data not available. 
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and N being the distance from the wall to the maximum in the local 
shear-stress profile. 

Thus using equations (3), (10) and (11) and scaling the turbulence 
structure on and A produces the result that 

P - pU s U m (12) 

or that 

AV/2 /«*u.V /3 

P ' X M \ h ) ' T M\N U m / (13) 

~~2 2 

Table 1 shows that p /t m values from experiments are in approximate 

agreement with Panton and Linebarger's calculations. It appears that 
~~2 2 

p /tjy| decreases with increasing H, although there are insufficient 

data to be sure of this and the a/N dependence given in equation (13). 

“ 2 2 

When t m /t w < 1.5, we would expect p /t m to approach 9, which is about 

the zero pressure gradient value. Hahn's (1976) data support this 

observation. The effect of transducer size on these data is not clear 
~2 2 

since p /t^ does not consistently decrease with dU^./ v . 

Figure 2 shows spectra corrected for transducer size, that 
correspond to the strong adverse pressure gradient data of Table 1 
when normalized on The data of Schloemer, Bradshaw, Burton and 
Lim for H < 2 and the calculations of Panton and Linebarger for a = 1 
agree within a few dB for w6*/U oo <5. All results seem to agree at 
<1)6*/^ =1, At higher values, Lims data are much higher than 

that of the others and Panton and Linebarger's calculations appear to 
be higher. As shown in Table 1, the u>’^ region becomes narrower at 
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' ) Figure 2: Spectra from some strong adverse pressure gradient 

turbulent boundary layers. Straight dashed line: 

(iT** variation. 





12 


higher H values until it is only a point near w6*/U s 0.5. A 

oo 

_3 

u> region at higher frequencies is observed for the higher H values. 

The a = 2 calculations of Panton and Linebarger agree fairly 

well with Burton's spectral data for H = 2.23. These calculations 

suggest that the ratios of streamwise length scale to length scales in 

other directions are greater than one for strongly retarded boundary 

layers. The spatial correlation data of Schubauer and Klebauoff (1951) 

for a strongly retarded separating turbulent boundary layer indicate 

that the streamwise integral length scale was about 2.5 times the normal - 

to-wall integral length scale. The calculations of Panton and Linebarger 
~2 p 

suggest that p /t m - a. The spectrum above = 1 does not appear 

to be influenced by a, either in calculations or in the data. Thus, the 
low frequency part of the spectrum for strongly- retarded adverse pressure 
gradient boundary layers appears to be large because a is much greater 
than one. 

Panton and Linebarger proposed that the convective velocity ll c (k) 

for the overlap region or k~^ part of the spectrum for an attached 

< 

boundary layer can be described by 


1) i u 

r ■ f «"/ k / + 5 -° 


(14) 


This expression is the same as the semi -logarithmic mean velocity profile 

equation, with U c = U at y = 1/k. They showed that this expression agrees 

very well with the data of Bradshaw (1967) and Wills (1970) for 6 < k6 < 40 

or the k~^ part of the spectrum. This expression also fits the measured 

U /U data of Schloemer within 10% for adverse, zero, and favorable 
c 00 
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minima. Panton and Linebarger suggested an approximate wake function 
expression for this low frequency range. 

In addition to these types of data, some measurements in zero 
pressure gradient flows have been obtained that lead to a better under- 
standing of the structural interrelationships between the velocity and 
surface pressure fluctuations. Thomas and Bull (1983) conditionally- 
sampled wall-pressure fluctuations on the basis of the high-frequency 
activity of the pressure fluctuations themselves, the high-frequency 
activity of the streamwise velocity fluctuations in the vicinity of the 
wall, and the excursions in velocity in the vicinity of the wall. This 
led to the identification of a characteristic wall-pressure fluctuation 
pattern which is associated with the burst-sweep cycle of events in the 
wall region. The pattern has the form of an overpressure over a stream- 
wise extent of about 1 '5-2-06* with a region of underpressure and a 
pressure minimum to either side of it, the distance between pressure minima 
being about 3-0-3-56*. This pattern is convected at a velocity 0-67 times 
the freestream velocity. Its phase relationship with velocity fluctuations 
close to the wall and the wall shear-stress fluctuations during the burst- 
sweep cycle were established. It appears to be produced by the inclined 
shear layer which forms the upstream surface of the large organized 
structures in the layer, and calculated pressure patterns support this 
conclusion. 

The phase relationships indicate that fluid involved in the bursting 
process is subjected to a favorable streamwise pressure gradient by the 
characteristic wal 1 -pressure pattern at the time that the lift-up of 
low-speed streaks in the wall region begins. In addition, order-of- 
magnitude estimates suggest that the adverse pressure gradients associated 
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with the characteristic pressure pattern, even if their phasing with 
streak lift-up were appropriate, would be insufficient to initiate 
the lift-up. It was concluded that the streamwise pressure gradients 
associated with the pressure patterns do not play an active role in the 
dynamics of the wall flow and are not the direct cause of the bursting 
process. Unfortunately, no studies of this type have been performed in 
the presence of strong adverse pressure gradients or separation. 

Mabey (1972) presented a rough correlation of rms pressure fluctu- 
ations and spectra for step-induced separation and reattachment flows. 
The length of the detached flow zone S is a good normalizing length 
that causes the shapes of the plots to be similar. In p'/q vs. X/S 
plots, p> increases rapidly after detachment, reaching an order of 
magnitude greater maximum just upstream of reattachment. At this maxi- 
mum, p'/q^ varied from 0.04 to 0.1, depending on the type of detached 
flow. In nF(n) vs. nS/U^ spectral plots, the peak frequency is located 
just below a nS/U^ of one. Since no data for the Reynolds shear stress 
-puv profiles are available, it is not possible to correlate further 
these data. 

Kiya et a 1. (1982) showed similar results for a forward-facing step 

flow: p'/q^ reached a maximum of 0.14 at X/L = 1; nF(n) reaches a 

_2 

maximum at nS/U^srO.? and decays proportional to n at higher fre- 
quencies. Near reattachment, p'/t^ = 10, which is an order of magnitude 
greater than the values given in Table 1 for attached flows. 


III. DESCRIPTION OF THE MIND TUNNEL AND THE TEST FLOW 


A. Wind Tunnel 

This facility is same one used in earlier work on steady and un- 
steady separating turbulent boundary layers (Simpson et al . , 1981 a, b, 
c, 1983 a, b). The mainstream flow of the blown open-circuit wind 
tunnel is introduced into the test section after first passing through a 
filter, blower, flow dampers, a section of honeycomb to remove the mean 
swirl of the flow, seven screens to remove much of the turbulence in- 
tensity, and finally through a two-dimensional 4:1 contraction-ratio nozzle 
to further reduce the longitudinal turbulence intensity while accelerating 
the flow to test speed. 

Figure 3 is a side-view schematic of the 8 m long, 0*91 m wide test 
section of the wind tunnel. The upper wall is adjustable such that the 
free-stream velocity or pressure gradient can be adjusted. The sidewalls 
are made of float plate glass to prevent laser signal dispersion, while 
the upper wall is made of Plexiglas. 

The active boundary- layer control system, which is described by 
Simpson, Chew and Shivaprasad (1980), was installed on the non-test 
walls of the test section to inhibit undesirable flow three-dimension- 
ality and to prevent separation. Because the static pressure in the 
test section is time-varying in unsteady experiments, no passive 
boundary- layer control can be used that depends on a steady test section 
pressure higher than the pressure outside the tunnel. Highly two- 
dimensional wall jets of high-velocity air are introduced at the beginning 
of each of the eight-feet long sections. At the latter two streamwise 
locations the oncoming boundary layer is partially removed by a highly 
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Figure 3: Sideview schematic diagram of the test section with the steady 
free-stream separating turbulent boundary layer (Simpson et al. 
1981a) on the bottom wall. The major division on the scales are 
10 in. Note the baffle plate upstream from the blunt leading 
edge on the bottom test wall and side- and upper-wall jet 
boundary-layer controls. 



Figure 


4: 


Free-stream velocity and pressure gradient distribution along 
the tunnel center-line. 


C p .2(P-P 1 )/pU_ 1 .l-(U.-U. j )i 


U ,=49-4 ft s' 1 . 
^1 




two-dimensional suction system. 

The inviscid core flow is uniform within 0-05% in the spanwise 
direction and within 1% in the vertical direction. The test-wall 
boundary layer is tripped by the blunt leading edge of the plywood 
floor, the height of the step from the wind-tunnel contraction to the 
test wall being 0-63 cm. Smoke can be introduced uniformly into the 
boundary layer just upstream of this trip for use with the laser 
anemometer. 

B. Summary of the Nature of This Steady Free-stream Separating 
Turbulent Boundary Layer 

All data presented here in the earlier works on this flow were 

1 0 

obtained at atmospheric pressure and 25 + j C conditions. Figure 4 
shows the free-stream velocity distributions obtained along the tunnel 
center-line using the single-wire probe. This distribution was repeata- 
ble within 2-9% over the duration of these experiments, which is only a 
little greater than the uncertainty in measuring the mean velocity with 
a hot-wire anemometer 2-4%). Figure 4 also shows the non-dimensional 
pressure gradient dC /dx along the centerline of the test wall. Hence 

r 

C s 2(P-P.)/pU 2 . =1 - (U /U .) 2 , where i denotes the free-stream 

P ! OO 1 OO OO ] 

entrance conditions at a distance x of 7.6 cm. A five-point local least- 
squares curve fit of Cp data was used at each streamwise location to 
determine this derivative. Just downstream of the location of the 
second wall-jet boundary- layer control unit (2.5 m), the slope of the 
static pressure gradient changes sign. Near 3.7 m the pressure gradient 
drops to an approximately constant value downstream. As discussed 
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in more detail by Simpson et al . (1981 a), the mean boundary- layer 
flow appears to be two-dimensional in the center third of the channel 
width upstream of 4.3 m. 

For steady free-stream mean- two-dimensional separating turbulent 
boundary layers, a set of quantitative definitions on the detachment 
state near the wall has been proposed (Simpson 1981): incipient 

detachment (ID) occurs with instantaneous backflow 1% of the time; 
intermittent transitory detachment (ITD) occurs with instantaneous 
backflow 20% of the time; transitory detachment (TD) occurs with 
instantaneous backflow 50% of the time; and detachment (D) occurs 
where the time-averaged wall shearing stress 7 = 0. Thus the fraction 

W 

of time with forward flow Yp^ is a descriptive parameter for identify- 
ing these stages and should be documented in all separated-flow experi- 
ments. 

Figure 3 shows a qualitative sketch of the steady free-stream 

bottom-wall turbulent shear flow studied with a laser anemometer at 

SMU and the locations of ID, ITD and D when determined 1 mm from the 

wall. The mean flow upstream of ID obeys the 'law of the wall' and the 

'law of the wake' as long as the maximum shearing stress - p'uv is 

less than l*5x . The qualitative turbulence structure is not markedly 
w 

di f ferent from the zero-pressure-gradient case. The 'bursting' fre- 
quency n of the most-energetic eddies near the wall is correlated by 
U,„/nc$ = 10, where U,,, is the mean velocity outside the boundary layer and 
■5 is the boundary-layer thickness. 

When - p77 > 1*5t , the Perry and Schofield (1973) mean- 

fTldX W 

velocity-profile correlation and the law of the wall apply upstream o 
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ITD. Up to one-third of the turbulence energy production in the 
outer region is due to normal -stress effects, which modify the rela- 
tions between dissipation rate, turbulence energy and turbulent 
shearing stress that are observed farther upstream. The spanwise 
integral lengthscale of the turbulence increases with 6 , and the 
bursting frequency n continues to be about equal to IU/106. Pressure- 
gradient relaxation begins near ITD and continues until D. 

Downstream of detachment, the mean backflow profile scales on the 
maximum negative mean velocity and its distance N 1 from the wall. A 

U + vs. y + law of the wall is not consistent with this result since U N 

1 /2 

and N increase with streamwise distance whilev/(t/p) / varies with 

W 

(t/p) 1 ^ 2 . High turbulence levels exist in the backflow, with u- and 
v-fluctuations of the same order as |U|. -uv/u'v' becomes lower with 
increasing backflow, and is about 25 % lower in the outer region than 
for the upstream attached flow. Mixing-length and eddy-viscosity models 
are adequate upstream of detachment and in the outer region, but are 
physically meaningless in the backflow. Yp U never reaches zero, indi- 
cating that there is no location with backflow all of the time. Normal 
and shear stresses turbulence energy production in the outer region 
supply turbulence energy to the backflow by turbulence diffusion where 
it is dissipated. Negligible turbulence-energy production occurs in 
the backflow. 

This turbulence-energy diffusion and the small mean backflow are 
supplied intermittently by large-scale structures as they pass through 
the detached flow as suggested by Figure 3. The backflow does not 
come from far downstream. The frequency of passage n of these large-scale 



structures also varies as U^/6 and is about an order of magnitude 
smaller than the frequency far upstream of detachment. Reynolds 
shearing stresses in the backflow must be modelled by relating them 
to the turbulence structure and not to local mean-velocity gradients. 

The mean-velocity profiles in the backflow are a result of time-averag- 
ing of the large fluctuations and are not related to the cause of the 
turbulence. 

IV. INSTRUMENTATION AND EXPERIMENTAL TECHNIQUES 

A. Hot-wire Anemometers and Probes 

Miller-type (1976) integrated circuit hot-wire anemometers and 
linearizers, as modified by Simpson, Heizer, and Nasburg (1979) were 
used. A standard TSI model 1244-T1.5 dual hot-wire probe was used 
for the stream-wise velocity fluctuation wavespeed and spectral 
measurements with a 5mm spacing between the parallel sensors. The fre- 
quency response for each sensor was flat up to 7.5 KHz for an overheat 
ratio of 1.7. The upstream sensor was used to obtain all spectral data. 

Calibrations were made in a TSI Model 1127 calibrator and in the 
core flow of the throat of the wind tunnel test section at various speeds. 
There was no detectable drift of the anemometer, although the function- 
module-type linearizers had a small amount of d.c. drift. Each linearized 
calibration had a low level of dispersion from a straight line, with the 
product moment correlation coefficient in excess of 0.9999. 

B. Microphones and Calibrations 

Two Sennheiser Model MKH-110 Instrumentation Microphones were each 
mounted in identical plexiglas housings with pinholes for use in these 
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measurements. According to the manufacturer' s specifications, these 
low frequency RF condensor microphones have a nominal sensitivity of 
20 mv/Pa, nearly flat frequency response (+ 3dB) between 1 Hz and 6 
KHz, an overload level of 120 dB SPL, and signal-to-noise ratio of 63 
dB + 3dB. The frequency response peaks at about 10 KHz and drops 
below the low frequency level of 16 KHz. Since the microphone can 
only be used with equal static pressures on both sides of the diaphragm, 
the back of the pressure sensititive element was sealed from the solid- 
state electronics portion. The sliding fit between the microphone and 
the plexiglas housing acted as a low-pass filter and permitted the volume 
behind the diaphragm to reach the test static pressure through vent holes 
on the side of the microphone. 

The 2.21 cm diameter cylindrical plexiglas housings each had 0.074 
cm diameter by 0.025 cm long pinholes in the center of the flat end. The 
1.55 cm diameter end of a microphone fits snugly against the housing end, 
with the 1.35 cm diameter diaphragm 0.10 cm away from the pinhole. 

Each microphone and plexiglas housing combination was calibrated 
with a Gen rad Model 1956 Sound-Level Calibrator between 125 Hz and 4 
KHz and 74 to 114 dB SPL. A calibrator consisting of a sealed volume 
with an oscillating piston was built and used for calibration between 
0.5 Hz and 30 Hz. Using the perfect gas law the adiabatic process equa- 
tions, the pressure oscillation in the sealed volume could be calculated 
at each test frequency from the volume oscillation produced by the piston. 
At low frequencies the measured sensitivity was approximately 22 mv/Pa. 

The following equation approximately describes the frequency response of 



these microphones in their housings 


P ”1 /2 

— = (1 - 0.437871 f 2 - 0.0800 f 4 + 0.0493292 f 6 ) 
p rms i 

where the frequency f is given in KHz. This equation indicates that 
the microphone sensitivity is reduced above 2 KHz because of the housing. 

For the coherence and celerity measurements when the two microphones 
are close together, two housings were joined adjacently. Several above- 
mentioned pinholes were located on each housing end, permitting sensing 
hole spacings of 1.40 to 3.02 cm. Pinholes not in use were covered by 
thin cellophane tape. 

C. Pressure Fluctuation Measurement Techniques 

After initial use of these microphones, it was apparent that minute 
wind tunnel vibrations and acoustic disturbances make large contributions 
to the microphone signals. To eliminate the vibration effects, the 
microphones and housings were supported directly from the concrete floor- 
ing, with no direct firm contact with the adjacent test wall. Strips of 
0.05mm thick celephane tape were used to make the test surface flat and 
continuous between the 3.5cm diameter holes in the test wall and the 
plexiglas housings. 

In this wind tunnel, there are streamwise acoustic disturbances that 
propagate downstream from the blower-muffler-plenum-honeycomb-screen 
sections into the contraction. The contraction and test section act as 
wave guides for these disturbances, so that at any streamwise location 
the streamwise acoustic waves are the same across the test section at any 
instant in time. The turbulent-flow-produced spectrum is the same across 
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the test section because the mean flow and mean square turbulence struc- 
ture are two-dimensional across the center of the flow (Simpson et al . , 
1980). The acoustic and turbulent signals are uncorrelated, since we 
can see from equations (2) and (3) that the turbulence-produced signal 
is due to the local ly-produced velocity field. These observations 
permit us to use two microphones to decompose the surface pressure 
fluctuation signals into the propagated acoustic part and the turbulent- 
flow-generated portion. 

The two microphones and housings are spaced far apart spanwise 
across the test section so that the turbulent signals are uncorrelated. 
The minimum distance between sensors to produce a zero cross-correlation 
is about 1/2 6, where <5 is the shear layer thickness (Simpson et al . , 
1977). Physically this means that individual large-scaled motions are 
no more than about 6 in spanwise extent and are unrelated to one another. 
In practice, for spectral measurements, aZ was 15 cm far upstream of 
detachment and 30 cm in the detached zone. 

At any spectral frequency n, the time-dependent signals detected by 
these two microphones can be written as 

Pin = Pi An + PlTn 05 a) 

P2n = P2An + p 2Tn 05 b) 

where the subscripts A and T denote the acoustic and turbulent portions. 
If we subtract p£ n from p-j n and obtain the mean square value of the 
result, then 



( p ln " p 2n^ 


2 


(16) 
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is the turbulent spectral contribution at frequency n because 


T "T 

p lTn = p 2Tn 


(mean 2-D flow) 


(17 a) 


p lAn p lTn = p 2An p 2Tn~ p lAn p 2Tn = p 2An p lTn = 0 * 17 

(uncorrelated acoustic and turbulent contributions) 


p lTn p 2Tn = ® (uncorrelated turbulent contributions) (17 c) 

p lAn = p 2An (same acoustic signals) (17 d) 

Under these conditions the acoustic contribution at frequency n is 
given by 


_ 2 

P 1 An 


(p ln + p 2n r 


p lTn 


(18) 


4 2 

The two microphones are located equi -distant about the tunnel 
centerline. Thus, equation (16) provides the proper turbulent spectrum 
for frequencies below c/w, where c is the speed of sound and w is the 
contraction and test section width (370 Hz). In other words, longitudinal, 
vertical and spanwise acoustic contributions that are the same at the 
two microphones are eliminated. However, spanwise acoustic contributions 
near frequency c/w and higher harmonics are added because of the anti- 
symmetry of that mode. For those frequencies the turbulent contribution 
must be obtained by 


(p 1n + P 2 „>‘ 


(19) 


HTn 2 

An experimental uncertainty analysis shows a negligible uncertainty in 
the resulting turbulent spectrum obtained in this manner when the sensi- 
tivities of the two microphones are matched, as discussed in section 



IV. D below. 

To determine the convective wave speed of the turbulent contributions, 
the two microphones were spaced a small distance aX apart in the streamwise 
direction. The wave speed or celerity at frequency n is given by 


U 


cn 


_ 2twaX 


(20) 


where 

tan ^n = V R n ’ ( 21 ) 

Y 2 (AX,n) = R n 2 + I n 2 (22) 

Because the acoustic contributions at two different streamwise locations 

are coherent but time delayed, they can be accounted for using the 

measured acoustic spectra. In practice these contributions were found 

2 

to be a very small part of y for frequencies where the coherence y 
produces meaningful and relatively certain results (Stegan and Van Atta, 
1970). 

D. Signal Processing 

For the pressure spectra and celerity data, the signals from each of 
the microphones were input to a TSI Model 101 5C Correlator where the 
effective volts/Pa sensitivities were equalized very closely with slightly 
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different amplification ratios. Thus, with closely equal effective sensi- 
tivities, it is possible to use the sum or difference of these signals to 
determine pressure spectral and celerity information, as discussed in 
Section IV. C above. Velocity spectra and celerities were obtained in the 
same manner by equalizing the effective volts/mps sensitivities of the 
two hot-wire anemometers. 

The sum or difference of these two instantaneous signals was input to 
a single-channel Princeton Applied Research (PAR) Model 4512 Fast-Fourier- 
Transform Spectrum Analyzer, The frequency resolution is 1/512 of the 
selectable frequency range of this 512 bin digital unit. Data records 
were at least 2 minutes long, which was experimentally verified to be 
sufficiently long for closely repeatable spectral results. A CRT display 
with a cursor permitted on-line viewing of spectral data. 

The logarithm of the content of each bin of a spectrum was transferred 
from the PAR to a HP 9825T digital computer for storage on tape and further 
processing. Use of the logarithms of the data preserved a 1 bit un- 
certainty of the spectral results for the final processing. Mean square 
fluctuation values were computed by summing the square of the anti -logarithm 
values for each bin, taking into account the proper transfer function and 
calibration values. 

V. EXPERIMENTAL RESULTS 

A. Streamwise Velocity Fluctuation Spectra 

Figures 5-12 show spectra for the streamwise velocity fluctuations 
upstream of detachment in the region of strong adverse pressure gradients. 
Simpson et al . (1980) report mean square velocity fluctuations for this 
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Figure 7: Velocity spectra for the downstream x location of 2.222 m. 

(The parameter 2tt6*/IL=0.001 7 sec., also note the displaced 
ordinates) 
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Figure 8: Velocity spectra for the downstream x location of 2.222 m. 

(The parameter 2ir6*/U»=0.0017 sec., also note the displaced 
ordinates) 


dB=10- 1og|n»F(n)| 
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Figure 11: Velocity spectra for the downstream x location of 3.009 m. 

(The parameter 2 Tr6*/U»= 0.0089 sec., also note the displaced 
ordinates) 
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flow. Here we present the spectra in the form of 10 log-|Q/nF(n)/ vs. 
log-|Q/n/ since 

oo 

f CnF (n ) ] d(tn/n/) = 1 (23) 

o 

by definition of F(n). Thus the peaks or plateaus of nF(n) plots show 
the most energetic spectral frequencies. 

Near the wall in the semi -logarithmic mean velocity profile region, 
nF(n) is nearly constant (n -0 * 1 ) over a relatively wide frequency range. 

A definite F(n) ~ n”^ inertial subrange spectral region is seen at 
higher frequencies as expected. As detachment is approached, smaller 
F(n) ' n~^ regions are observed. No specific frequencies with large 
contributions are observed, so the wind tunnel acoustic fluctuations are 
definitely not interacting with the turbulence to produce velocity spectrum 
contributions. At progressive downstream locations, the frequency range 
for the weak contributions decreases, as Simpson et al . (1977) observed 

for a slightly different separating flow. Near and downstream detachment, 
the inertial subrange begins at progressively higher frequencies at higher 
y/6 for a given X, as observed in Figures 12-16. The U^/nfi values for the 
peak frequency of each nF(n) distribution in the middle and outer boundary 
layer regions are about 10 +3, as observed by Simpson et al . (1981 b) 
upstream and downstream of detachment from laser anemometer data. 

B. Streamwise Velocity Fluctuation Celerities 

Figures 17 and 18 show some velocity fluctuation celerity results 
calculated using equations (20-22) for the data obtained from the parallel 
hot-wires that are a distance AX = 5 mm apart. The wave celerity U cn is a 
function of spatial position and frequency, but not the streamwise spacing 
between the two sensors (Favre et al . , 1967). As observed by 
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Figure 15: Velocity spectra for the downstream x location of 3.970 m. 

(The parameter 2ir«*/U»*0.0734 sec., also note the displaced 
ordinates) 
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Figure 16: Velocity spectra for the downstream x location of 4.340 m. 

(The parameter 2tt6*/Uoo = 0.0996 sec., also note the displaced 
ordinates) 
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Fit of the model equation for each y/6 position 


Figure 17: Velocity celerity (U c ) versus frequency (Hz) for the 
downstream x location of 3.009 m. 

(Note the displaced ordinates) 



U c (m/sec) 






r.X3im 


0.4922 
0.5778 
0.6526 
0.7382 
0.8132 
0.8988 

(From top to bottom) 







\um 


jMu 




0.4 i.o 


10.0 

i (Hz) 


100.0 


Fit of the model equation for each y/s positii 


Figure 18: Velocity celerity (U c ) versus frequency (Hz) for the 
downstream x location of 3.970 m. 

(Note the displaced ordinates) 
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Strickland and Simpson (1973) for their separating turbulent boundary 
layer, the apparent wavespeed increases with increasing frequency. 
Strickland and Simpson obtained U cn = aX/t values at each frequency, 
where t is the time delay for the maximum cross-correlation of the two 
narrow bandpass-filtered siqnals. The data are shown for coherence 
values above 0.36, which Heidrick et. al. (1971) found to be the minimum 
value that reasonably reliable results could be obtained. Figure 19 
shows the values of 2irnAX/U for each y/s at which the coherence is 0.36. 
For increasing y/6, ZirnAX/U is lower for the same coherence. 

At high frequencies we expect that the small est-scale motions move 
on the average with a wavespeed equal to the local mean velocity. If 
the large structures move at a lower velocity, then the defect with the 
local mean velocity should be proportional to the wavelength. In other 
words. 


U - U U 
cn _ cn 

U " 2nnL 


(24) 


where Lis some length scale describing the flow at that spatial position. 
This equation can be rearranged to produce 

^cn 2irnL/U 

U " 2itnL , , (25) 

U 1 

which is shown in Figure 20 and will be referred to as the 'model' equation. 

Figures 17 and 18 show crude curve fits of the model equation to the 
lower frequencies. Since U at each position is known, the only adjustable 
parameter in equation (25) is L. Figures 21 and 22 show the experimental 
values of L/s and L/s* determined at all available spatial locations. 
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+ x=2.700 m X x=2.854 m 

* x=3.009 m □ x=3.228 m 

O x=3.527 m A x=3.679 m 

Y x=3,970 in Z x=4.340 m 


Figure 19: y/6 versus 2TrnAx/U where R n =0.6 
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Figure 20: Model equation versus 2irnL/U 
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Legend 
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Figure 21: L/6 versus y/6 

(Note the log-linear abscissa) 





These L values are about + 30% uncertain due to the subjective nature 

of fitting equation (25) to the mid-range of frequencies. 

L decreases with increasing streamwise position, reflecting the 

strong deceleration of the large scale structures. At a given stream- 

wise position, L increases to a maximum, which lies increasingly farther 

from the wall at downstream locations. Downstream of detachment (3.52 m) , 

no U measurements near the wall could be obtained because of the inter- 
cn 

mittently fluctuating flow direction. Low values of U cn are observed at 
y/6 = 0.5, the nearest to wall location without intermittent backflow. 
Consequently, L is very low in these nearest to the wall data. Figure 23 
shows that L/aX is not a detectable function of 2nnAX/U when the 
coherence is 0.36, which indicates that the level of the coherence for 
these data does not strongly influence the results for L. 

C. Surface Pressure Fluctuation Spectra 

Figure 24 shows surface pressure fluctuation spectra nF(n) for the 

adverse pressure gradient region upstream of incipient detachment (ID) 

at 3.11 m where there is no intermittent backflow near the wall. As 

documented by Simpson et al . (1981 a) for this flow, the wall shear 

stress decreases along the flow while the maximum shear stress decreases 

~~2 

slightly. As shown in Table 2, T max / T w and p increase along the flow. 

Figure 25 shows surface pressure fluctuation spectra downstream of 
the beginning of intermittent backflow. The maximum shear stress increases 
while the mean wall shear stress decreases to zero at detachment at 3.52 m, 
and is negative downstream. Figure 26 shows that the bandwidth of the 
plateau in the nF(n) curves decreases with streamwise distance until it is 
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Figure 23: L/ax versus 2TtnAx/U 
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Figure 24: Surface pressure spectra upstream of incipient detachment 
(Note the displaced ordinates) 
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no longer detectable downstream of incipient detachment. 

The results in Figures 24 and 25 were obtained by subtracting the 
signals from the two microphones that have equal effective signal 
sensitivities, i.e., equation (16) and the technique described in section 
IV. C were used. For frequencies below 15Hz, F(n) is about constant up- 
stream of detachment. For higher frequencies less than 150 Hz, nF(n) is 
about constant upstream of incipient detachment. 

Between 150 Hz and 1000 Hz two effects were observed. The vertical 
and spanwise acoustical waves that are anti -symmetric about the tunnel 
centerline are additive when the signals are subtracted and equation (16) 
is used. Consequently, equation (19) was used to obtain the turbulent 
contributions at those frequencies. The straight lines shown on Figure 
24 show the approximate corrected spectrum for that frequency range. 

The second effect is a hump in each spectrum. At first glance this 
is apparently the pinhole-related amplification of a certain frequency 
range described by Bull and Thomas and discussed in Section II. This 
effect does not cancel when subtracting the signal from two equally sensi- 
tive microphones when using equation (16). Bull and Thomas only examined 

zero-pressure-gradient cases and concluded that ^ > 0.1 for this effect 

IT 

T 

2 

to be important. On this criterion, wv/ll^ < 0.1 for the present data 
under 1 KHz upstream of incipient detachment, making this pinhole effect 
appear to be unrelated to the humps in the nF(n) spectra near 500 Hz. 

If we re-examine the Bull and Thomas data, we see that the wave- 
lengths that were amplified in their experiments are equal to the pinhole 
diameter when 1 < u>v/U T < 2. In other words, when nd/U c =;l, spurious 
amplification of pressure fluctuations occurred for the Bull and Thomas 
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pinhole data. The present data produce values of nd/U < 0.1 for 
the observed humps. Thus, it is reasonable to conclude that the Bull 
and Thomas pinhole effect is unimportant for these data and that the 
observed humps are produced by the turbulent flow. If one compares 
these spectra with the velocity spectra in Figures 5-11, it appears 
that the low frequency end of each hump coincides closely with the low 
frequency end of the n region for the nF(n) velocity spectra. 

Figures 27-29 show the spectral data for the attached flow in 
non-dimensional coordinates. The Corcos (1963) sensor resolution cor- 
rection was applied to these data and amounted to the order of 1 dB at 
the higher frequencies. The freestream dynamic pressure has less direct 
influence on the pressure spectral than the local streamwise maximum 
shearing stress t^. Consequently, Figure 28 shows a tighter correlation 
of the ordinate values than shown in Figure 27. In both of these figures, 
the higher frequency portions of the spectra (u^i/U^ > 1) do not correlate 
well. Figure 29 shows a much better correlation of spectra at all 
locations upstream of detachment when t m is used to non-dimensionalize 
the frequency as well as |(io). The variable wv/U M for strong adverse 
pressure gradient cases is analogous to ojv/U^ , which correlates data 
at the higher spectral frequencies for zero pressure gradient flows. 

In fact, for zero and favorable pressure gradients * U ; As shown 
in Table 2, there is not much variation of l^/U^ along the attached 

flow, SO idv/U^ ~ (wS-j/Uj (yg-) . 

1 °» 

Figure 28 shows that for a)6-|/U oo < 1, the spectra varies about like 
of 0 - 7 , which is observed in the data of Brooks and Hodgson (1981) for 
a NACA 0012 airfoil and was observed by McGrath and Simpson (1985) for 
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zero and favorable pressure gradient flows with Re„ < 5000. Upstream 

U 

of the strong adverse pressure gradient region (X < 2.54 m and t^/t < 1.5), 

-8 

the higher frequency spectra vary like u> . McGrath and Simpson (1985) 

-5 5 

observed a w ' variation for the highest frequency components. In 

-3 

the strong adverse pressure gradient region, <f> (u>) ~ w at the higher 
frequencies, which is not shown in Figure 24 but was observed from data 
obtained using the 2 KHz range of the FFT that are included in Figures 
27-29. 

-2 

Figure 25 clearly shows the n dependence of nF(n) at the higher 

_? 

frequencies of the detaching flow. This n behavior has been observed 
in the data of Bull and Thomas for a zero-pressure-gradient flow, Kiya 
et al . for a separated forward facing step flow, and Bradshaw, Burton, 
and Schloemer for adverse pressure gradient boundary layers. When we 
compare Figure 28 and Table 2 with Figure 2, we see that $>((u) has 
approximately the same shape and magnitude as for the earlier attached 
flow investigations. 

Figure 25 shows that the frequency range where F(n) is constant 

decreases continuously in the downstream direction. At the lowest 

frequencies downstream of detachment, F(n) varies with n to a power 

greater than one. Figure 30 shows that the peak frequency of nF(n) and 
_2 

the n high frequency range of the spectra in Figure 25 can be correlated 

using the non-dimensional frequency a>6*/U oo . The peak occurs at 

_2 

a)6*/U = 0.8 and the n range extends to about 6. Figures 12-16 show 

oo 

that the peak frequencies of nF(n) for the velocity spectra near the 
maximum shearing stress location (y/6 ~ 0.5) are approximately the same 
as peak frequencies shown in Figure 25. The streamwise velocity spectra 
in the intermittent backflow region obtained by Simpson et al . (1981 b) 
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Figure 27: Surface pressure spectra upstream of incipient 

detachment, normalized on local freestream q, 

6,, and U : # , 1.63 m; ▲, 1.89 m; T , 2.22 m 

I 00 

■ , 2.54 m; ♦ , 2.85 m. 






10 Log 



0.001 0.01 0.1 1.0 10.0 

tub 1 /Uoo 


Figure 28: Surface pressure spectra upstream of incipient 

detachment, normalized on local x M , 6., and U . 
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variation; broken line: <d variation. 
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Figure 29: Surface pressure spectra upstream of incipient 

detachment, normalized on t m - 



dB=10»log|n*F(n)| 
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Legend 

+ x=3. 112 m 

X x=3.276 m 
* x=3 . 421 m 
□ x=3.524 m 
O x=3.749 m 
A x=4,140 m 
Y x=4.419 m 


Figure 30: Nondimensional surface pressure fluctuation spectra 

downstream of incipient detachment. Broken line: u> 2 - 4 ; 
Dashed line: w 1 ; Solid line: u>~ 2 . 
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with a laser anemometer indicate that nF(n) is constant for 0.07 < 

JR ^ 

£j — <0.5. This frequency range is below the peak frequency in 

oo 

Figure 30. 

Figures 31 and 32 show the surface pressure spectra downstream of 
incipient detachment in non-dimensional coordinates. In Figure 31 
$(u >) is normalized on q; this plot masks the true variations at low 
frequencies that are observed in Figures 25, 30 and 32. Figure 32 
shows that normalization on x M produces a much tighter correlation for 
U)6 l/ U » > ^ • The Corcos (1963) microphone resolution correction was not 
applied to these data because large wavelength motions dominate the 
detached flow and require negligible correction. Furthermore, and more 
importantly, the instantaneous flow reversal contains small and large 
wavelength motions that as yet cannot be related to specific frequencies. 

Figure 33 shows the rms pressure values p' computed from these data 
and normalized on the reference inlet dynamic pressure, the local wall 
shear stress x , and the local maximum shear stress x M - The uncertainty 
of p' is about + 20%. Although p’ increases along the flow, p'/x^ increases 
to detachment and then decreases. At X=1 .63 m at the end of the favorable 
pressure gradient region, p'/x =2.3 which is in very good agreement with 

w 

values obtained by McGrath and Simpson (1985) for a favorable pressure 
gradient flow at the same Re 0 . 

D. Surface Pressure Fluctuation Celerities 

Figures 34 and 35 show the square root of the coherence for two 
microphones whose pinholes are separated aX in the streamwise direction. 
These results were obtained using equations (20-22). The values of y 






Figure 32: Surface pressure spectra downstream of incipient detachment, normalized 

on t m . Broken line: w ^ variation. 
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Figure 33: Streanwise distribution 
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are low, indicating a large decay in similarity of the pressure signals 
between microphones. Like the data of Brooks and Hodgson and of Hahn, 
y increases with frequency to a maximum in the frequency range where 
nF(n) is a maximum for the locations upstream of incipient detachment. 

At higher frequencies y drops to very low values. 

Corcos (1963) proposed that the cross spectrum in the lateral 
and longitudinal directions decay exponentially with the phase angle <j>. 

In terms of the square root of the coherence 

Y * exp (-K-|<j> - K 3 ^fj^,+ ) (26) 

Upstream of incipient detachment ^=0.12, as shown in Figure 34. At the 
lower values of <t> , agreement with equation (26) is poor, but Brooks and 
Hodgson show a similar behavior. This value of seems low for a strong 
adverse pressure gradient flow in view of Schloemer's (1967) result that 
K.| is greater for adverse pressure gradients than for zero pressure 
gradients. Zero pressure gradient values of 0.11 to 0.23 have been 
reported (Corcos, 1963; McGrath and Simpson, 1985; Brooks and Hodgson, 
1981 ). 

After some intermittent backflow begins downstream of incipient 
detachment, the value of K-j appears to drastically increase, as shown 
in Figures 34 and 35. A K^=0.75 approximately describes the data at 
X=3.28 m and 3.53 m; K-j - 1.0 at X=4.14 m. These values of are close 
to the values of near 0.75 reported for the lateral or spanwise decay 
of y for zero and adverse pressure gradient attached boundary layers 
(Corcos, 1963; McGrath and Simpson, 1985; Brooks and Hodgson, 1981; 
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Figure 34: Square root of streamwise coherence upstream of 

and in the detachment zone. 

Upper line: K 1 = 0.12 in equation (26) 

Lower line: K-| =0.75 in equation (26) 

O » 2.22 m, aX/6-j = 2.55; 0 » 2.85 ro> aX/ *S ^ = ( 

3.28 m: V , AX/6, = 0.28; a, aX/6, = 0.44; 
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Figure 35: Square root of streamwise coherence downstream 

of detachment. Top line, equation (26) with 
K-| = 0.75; Bottom line, equation (26) with = 

At 3.53 m: O, AX/6 ] = 0.177; □, aX/6 ] = 0.28 

0 , AX/6^ = 0.383. At 4.14 m: •, aX/6^ * 0.0 

■ , AX/6, = 0.122; ♦, AX/6, = 0.165. 
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Figure 36: Pressure fluctuation celerity: ▲ , 2.22 m; 

• , 2,85 m; ■, 3.28 nr, ♦, 3.53 m; 0, 

4.14 m. Solid lines: results from equation 

(14) for first three locations. Note displaced 
ordinates. 
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Blake, 1970; Burton, 1973). This suggests, as does the flowfield data, 
that the intermittent backflow destroys the streamwise coherence of 
pressure fluctuations. 

Figure 36 shows that upstream of incipient detachment, U c increases 
with increasing frequency till near u^i/U^ = 1 as observed by Brooks and 
Hodgson. U decreases at higher frequencies and agrees asymptotically 

V 

with equation (14) for the inner-outer overlap region. After the begin- 
ning of intermittent backflow, U c no longer agrees with equation (14) 
and very low values of U c are observed. The experimental uncertainty on 

U is about ± 0.05. 
c 

VI. DISCUSSION, SUMMARY AND CONCLUSIONS 

The streamwise velocity fluctuation spectra and celerities for this 
flow possesses expected features. There is no effect of wind tunnel 
acoustic noise on the velocity spectra. Near the wall upstream of 
detachment, nF(n) is nearly constant over a frequency range while a 
definite F(n) - n~ 5 ^ 3 inertial subrange region is observed at higher 
frequencies. The frequencies of the energy-containing motions decrease 
as detachment is approached, as observed by Simpson et al . (1977, 1981 b). 

As observed by Strickland and Simpson (1973), the celerity U cn 
increases with frequency till it approaches the local U asymptotically. 

The simple crude model equation (25) was observed to fit U cn along this 
flow. Unfortunately, the length scale L in this model equation is not 
known a priori. It was observed to decrease with streamwise position 
and to reach a maximum at a given streamwise position increasingly farther 
from the wall. Downstream of the beginning of intermittent backflow, no 
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U cn measurements near the wall can be made because of the ambiguity of 
the flow direction. 

At the end of the favorable pressure gradient region far upstream 

of incipient detachment, p7x w = 2.3, which is in very good agreement 

with values obtained by McGrath and Simpson (1985) for a favorable pressure 

— J 2 

gradient flow at the same Re 0 . The ratio p /x M is substantially larger 

for the attached strong adverse pressure gradient part of the present flow 

than in the investigations listed in Table 1. As discussed in Section II, 
~~2 2 

p /t^ is proportional to the ratio of streamwise length scale to length 

scales in other directions a. The integral length scale data of Chehroudi 

and Simpson (1983) show that a is about 2.5 for this flow, which is in 

agreement with the Schubauer and Klebanoff (1951) results for a in their 

~2 2 

Separating flow. Thus, the present results for p /t m upstream of 

~2 

detachment appear plausible. This apparent importance of a on p suggests 

i 

that space-time velocity correlations should be obtained in future experi- 

e 

ments to determine integral length scales. 

~2 

Downstream of detachment, p continues to increase, although p /q re f 
does not reach values of 0.04 to 0.1 that were observed by Mabey (1972) 
for step-induced separation and reattachment flows. The ratio P'/t^ 
increases to detachment and then decreases downstream. The level of 
P'/t m for these data seems reasonable since Kiya et al . (1982) observed 
values of p'/t^ s 10 for a forward-facing step flow. This decrease 
appears to be due to the fact that the pressure-fluctuation-producing 
motions move increasingly away from the wall downstream of detachment. 

Both the turbulence-mean shear interaction and the turbulence- 
turbulence interaction in the pressure fluctuation source term, equation 
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(2), are important for detached flows. Velocity fluctuations are 

as large as mean velocities in the backflow. Reynolds shear stresses 

and their gradients are large away from the wall. Thus the largest 

pressure fluctuations are not at the wall in a detached flow, but must 

be near the middle of the shear layer. Equation (3) indicates that 

p(x) on the wall must decrease if the source q moves away from the wall. 

Table 2 shows that the distance from the wall to the maximum shear 

location, N, increases rapidly downstream of detachment. 

In order to try to correlate this effect a modified version of 

equation (13) was used. First, N/6* is almost constant downstream of 

1/3 

detachment and (Uyil^)' does not vary enough to correlate p' / t^ . As 
shown in Table 2, i , 



(27) 


where is the distance of the maximum shear stress location from the 
wall at detachment. 

~2 

The spectra for p substantially agree with other investigations. 
Upstream of incipient detachment, u>4>(w) is nearly constant over a 
frequency range that decreases as detachment is approached. At lower 
frequencies $(u>) varies like at higher frequencies $(<u) varies 

_3 

with w . Downstream of the beginning of intermittent backflow, oi<t>(oj) 

_o 

has a peak near = 0.8 and a w range at higher frequencies. 

For lower frequencies, w<i>(u)) varies with u near the beginning of 

intermittent backflow; farther downstream this low frequency range 
2 4 

varies with u> ' . 
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The coherence of the pressure fluctuation producing motions 
remains high in the streamwise direction upstream of incipient detachment 
but drops drastically with the beginning of intermittent backflow. The 
streamwise coherence level looks much like that for the spanwise direction 
for attached flows. 

At low frequencies, both upstream and downstream of detachment, U cn 
celerity results are reasonable and agree qualitatively with other 
available data. Upstream of detachment at high frequencies the data 
agree asympotically with the overlap region celerity equation (14). 

The relatively large size of the two microphones used in this work 
prevented celerity measurements at closer spacings. As discussed in 
Section II, one should obtain U (k,oj) in future work by obtaining the 

v| 1 

Fourier transform of frequency-filtered spatial correlations, as done by 
Bradshaw (1967) and Wills (1970). Downstream of the beginning of inter- 
mittent backflow, the instantaneous wavespeed U cn can be both positive 
and negative for sufficiently high frequencies. Thus, the long-time- 
averaged U is lower than upstream for these frequencies. 

From the perspective of using pressure fluctuation data to calculate 
farfield noise, one should probably locate the effective pressure 
fluctuation sources along or near the locus of the maximum shear stress 
position. The pressure-fluctuation-producing motions are concentrated 
in this region. The celerities of these fluid motions are probably close 
to those for the pressure fluctuations in this region. 

This hypothesis is consistent with the idea of placing the effective 
pressure fluctuation sources on the flow surface for low pressure gradient 
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flows, since in these cases the maximum velocity gradients and shear 
stresses are at or very close to the wall. Since it is not possible 
to accurately measure pressure fluctuations within the turbulent 
separated flow, measurements of the farfield pressure fluctuation and 
the wall measure fluctuation should be used to estimate the effective 
pressure fluctuations at the maximum shearing stress location. 
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